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1  a  ABSTRACT  (M«mjn  200  urorda) 

A  reconf igurable  antenna  concept,  consisting  of  a  semiconductor  illuminated  by  an 
optical  or  electron  source,  was  studied.  With  appropriate  illumination,  mobile 
charge  carriers  excited  from  the  valence  band  to  the  conduction  band  can  form  a 
conducting  antenna  pattern  in  the  semiconductor  material.  Calculations  of  gain  and 
radiation  resistance  for  linear  dipole  and  "bow-tie"  layouts  were  completed  to 
estimate  source  power  requirements  and  practical  wavelength  limits  for  such  an 
antenna.  Significant  findings  were  that  source  power  for  writing  a  reconf igurable 
antenna  Increases  as  the  square  of  the  radiation  wavelength  ;  and  that 
semiconductor  material  permittivity  limits  antenna  efficiency  to  about  one-tenth  that 
of  a  comparable  free-standing  antenna.  Frequency  agility  in  the  antenna  may  be 
possible  using  an  electronically  controllable  liquid  crystal  display,  semiconductor 
laser  array,  or  electron  emitter  tips. 
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The  Reconfigurable  Antenna  and  Objective  of  Study 


The  reconfigurable  antenna  under  study  consists  of  a 
semiconductor  illuminated  by  an  optical  or  electron  source. 
The  illumination  is  non-uniform  over  the  semiconductor.  The 
illuminated  part  becomes  conducting  and  forms  the  antenna. 
The  conducting  pattern  can  be  controlled  electronically, 
making  the  system  agile.  With  light,  the  pattern  can  be 
illuminated  with  semiconductor  laser  arrays  whose  driving 
currents  are  controlled  electronically,  or  with  an 
incoherent  source  through  a  liquid  crystal  display  matrix 
whose  transmission  is  controlled  by  a  computer.  With 
electrons,  emitter  tips  can  be  controlled  with  electrical 
signals.  In  both  cases,  mobile  charge  carriers  are  created 
in  the  semiconductor  by  exciting  electrons  from  the  valence 
band  to  the  conduction  band.  The  excitation  energy  needed 
per  electron  is  the  bandgap  energy.  The  key  question  is 
whether  usable  antennas  can  be  written  with  practical  source 
(optical  or  electron)  powers,  and  over  what  frequency  range. 
This  study  endeavers  to  answer  this  question.  A  summary  of. 
the  findings  of  this  study  follows  immediately. 
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Summary  of  Findings  of  this  Study 

The  study  finds  that: 

(1)  The  source  (optical  or  electron)  power  needed  to  write  a 
reconf igurable  antenna  on  a  semiconductor  dncTssses  with  the 
sgu3T6  of  the  zsdistion  wavelength ,  thus  exacting  a  heavy 
price  on  low-frequency  operations.  For  silicon,  the  optical 
power  needed  to  write  an  antenna  one  wavelength  long  at  1 
GHz  is  about  300  watts,  and  electron  excitation  requires 
about  1  kW  for  the  same  antenna; 

(2)  Because  the  permittivities  of  semiconductors  are  usually 
10  times  that  of  vacuum,  an  antenna  on  semiconductors, 
reconf igurable  or  metallic,  is  about  10  times  less  efficient 
as  a  free-standing  antenna  of  the  same  electrical  structure; 

(3)  Because  of  current  decays  in  the  necessarily  resistive 
reconf igurable  antennas,  a  single  antenna  of  dimensions 
larger  than  a  wavelength  is  hard  to  realize.  For  high 
gains,  an  array  of  smaller  antennas  has  to  used; 


(4)  Wire-type  antennas  are  suitable  for  directional 
excitation  sources  like  lasers  and  micro-tip  emitters;  large 
aperture  antennas,  for  flood  sources  like  the  flash  lamp. 


(5)  For  reasons  stated  above,  practical  reconf igurable 
antennas  are  limited  to  elementary  structures  for 
applications  above  a  GHz  or  so,  if  the  antenna  systems  are 
to  be  inexpensive  and  simple. 
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Establishing  Criterion 


The  conductivity  of  patterns  written  on  a  semiconductor  by 
light  or  electrons  is  proportional  to  the  light  or  electron 
power  incident  on  the  semiconductor,  and  in  practice  cannot 
be  made  as  high  as  that  of  metal.  It  means  then  that  the 
reconf igurable  antenna  cannot  be  as  effective  as  the 
metallic  one,  but  can  only  approach  the  metallic  in  the 
limit  of  infinite  source  power.  A  criterion  must  be 
established  to  determine  the  feasibility  of  reconfiguration. 
There  are  two  obvious  parameters  to  determine  the 
performance  of  an  antenna:  the  gain  and  the  radiation 
efficiency.  Two  antennas  of  very  different  efficiency  can 
have  the  same  gain,  for  example,  the  very  inefficient  short 
Hertzian  dipole  has  almost  the  same  gain  at  all  directions 
as  the  much  more  efficient  half-wave  dipole.  Therefore  for 
this  study,  we  will  use  efficiency  as  measured  by  the 
radiation  resistance,  and  take  the  required  source  power  as 
the  power  needed  to  induce  sufficient  conductivity  so  that 
the  reconf igurable  antenna  radiates  half  as  effectively  as  a 
metallic  antenna  of  the  same  structure.  By  reciprocity,  the 
same  antenna  used  for  reception  will  be  half  as  effective  as 
the  metallic  reference.  Beyond  this  point  of  -3  dB 
efficiency,  calculations  reported  below  indicate  that  the 
radiation  effectiveness  increases  much  more  slowly  as  source 
power  increases. 
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Estimate  of  Required  Source  Power  and  Scaling 

An  order-of-magnitude  estimate  of  the  minimum  required 
optical  or  electron  power  for  writing  a  reconf igurable 
antenna  can  be  made  as  follows.  Most  radiating  systems  have 
characteristic  impedance  of  50  fl  to  100  Q.  For  microwave 
structures  on  dense  dielectrics/  characteristic  impedances 
are  reduced  by  the  refractive  index,  typically  3  to  4  in 
semiconductors,  therefore  the  resistance  of  the  antenna 
pattern  cannot  be  much  higher  than  about  20  O,  or  mismatch 
alone  will  render  the  antenna  ineffective.  The  minimum 
source  power  can  be  estimated  to  be  the  power  needed  to 
reduce  the  resistance  of  the  antenna  pattern  on 
semiconductor  to  20  Q.  Consider  electrons  or  light  of 
power  P  penetrate  to  depth  d  in  the  semiconductor  to  create 
an  antenna  strip  of  width  w  and  length  L,  current  flowing 
along  L  (Fig.  1) .  If  a  is  the  fraction  of  P  absorbed  by  the 
semiconductor  for  exciting  electrons  from  the  valence  to  the 
conduction  band  across  the  bandgap  Eg,  and  x  is  the  electron 

lifetime,  then  the  electron  density  n  can  be  calculated  by 
balancing  power  to  be: 

n  =  aPx/ (EgLdw)  ,  (1) 

whence  the  resistance  R  along  L  is: 

R  =  L2Eg/  (eiiaPx)  (2) 

where  is  the  electron  mobility  and  e  the  electron  charge. 
Any  antenna  of  sufficient  gain  and  directivity,  which  rougly 
scale  as  (size/^)2,  must  have  a  size  of  at  least  one 
wavelenght  X.  Equating  L  with  X.  in  Eq.  (2)  yields  the 
minimum  source  power  required: 

P  =  {e\iax)  (3) 
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The  significance  of  this  equation  is  that  the  required  power 
increases  with  the  square  of  the  wavelength.  The  wavelength 
X  is  the  wavelength  in  semiconductor,  typically  about  0.3  of 
that  in  vacuum. 


For  light  sources  in  the  visible  and  near  infrared,  the 
photon  energy  is  close  the  bandgap  energy  Eg,  so  a  ~  1 

except  for  narrow  bandgap  materials  like  germanium.  For 
low-energy  (<  few  keV)  electrons,  typically  3Eg  [1]  is 


needed  to  excite  one  electron  to  the  conduction  band,  and  a 
~  0.3.  Using  established  material  parameters,  the  minimum 
required  powers  are  listed  in  Table  1  for  silicon  and 
germanium.  (The  carrier  lifetime  of  silicon  depends  on  many 
factors  -  impurities,  defects,  processing  steps,  and  can 
vary  from  under  1  ps  to  1  ms.  The  more  typical  value  of  20 
|is  is  used  in  Table  1.  Also,  because  of  the  dielectric 
mismatch  between  vacuum  and  semiconductors,  about  30%  of 
incident  light  is  reflected.) 


Table  1 

Minimum  Source  Power  for  Reconfigurable  Antennas 


Optical  Source  (watts) 
Electron  Source  (watts) 


Silicon  Germanium 

0.04^2 

lOX^  0.06X2 

Wavelength  X  in  cm 


For  example,  for  1  GHz  in  silicon,  X  ~  10  cm,  the  minimum 
required  optical  power  is  300  Watts,  which  we  consider 
practical . 

The  estimates  given  here  are  validated  by  more  accurate 
calculations  shown  below  for  several  cases. 
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Calculations  on  Resistive  Antennas 
(i)  General  Discussions 

For  more  accurate  estimates  of  the  required  source  power, 
calculations  have  been  performed  on  two  representative 
antenna  structures  with  finite  electrical  resistivity.  The 
antennas  are:  (1)  linear  dipole;  and  (2)  "bow-tie."  The 
linear  dipole  is  the  simplest  useful  antenna,  out  of  which 
other  more  complicated  antennas  of  higher  directionality  and 
gain  can  be  constructed.  Directional  sources  like  lasers 
and  electron  beams  can  selectively  write  linear  dipoles 
efficiently.  The  bow-tie  is  a  large-aperture,  broad-band 
antenna  suitable  for  flood  sources  like  the  flash  lamp. 
Since  the  cross-sectional  area  of  the  conducting  path  of  the 
bow-tie  increases  as  the  distance  from  the  terminal 
increases,  the  resistance  decreases  as  the  current  flows 
towards  the  end  of  the  antenna,  which  to  some  degree 
counter-balances  the  effectively  reduced  antenna  length  due 
to  current  decay,  which  is  accomplished  without  non-uniform 
illumination . 

Throughout  the  calculations,  we  have  followed  the  standard 
simplification  of  ignoring  the  reaction  of  the  radiation 
field  on  the  current  distribution  [2] .  We  assume  an 
excitation  current  of  given  amplitude  and  frequency  feeding 
the  antenna,  the  current  propagates  along  the  resistive 
antenna  and  radiates.  The  total  radiated  power  divided  by 
the  current  squared  yields  the  radiation  resistance,  the 
measure  of  radiation  efficiency.  By  reciprocity,  the 
radiation  resistance  also  measures  the  reception  efficiency 
of  the  antenna.  The  radiation  impedance  of  the  resistive 
antenna  is  different  from  that  of  the  conventional  antenna. 
We  have  ignored  the  additional  problem  of  impedance 
matching,  since  it  is  only  incidental  to  the  basic  question 
of  power  source  requirement.  In  any  case,  the  difference  is 
about  a  factor  of  2  in  the  final  radiation  efficiency. 
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In  general,  the  more  exact  calculations  shown  below 
confirmed  the  1st  order  estimate  presented  in  the  last 
section. 

A  resistive  line  is  modeled  by  a  distributed  series 
resistance  per  unit  length  R'[2],  in  addition  to  the 
standard  series  inductance/length  L'  and  shunt 
capacitance/length  C  (Fig.  2).  All  time-dependence  is 
expressed  by  e^®^,  where  ©  is  the  angular  freguency  of  the 
radiation.  The  propagation  constant  y  along  the  resistive 
line  is  complex: 

Y  =  (joC'R'  -  ffl2L.ci)l/2 


=  jk(l  +  Xi)  ^/2exp(-0.5jtan~^xi)  (4) 

where  k  =  ©(L'C)^'^^  is  the  propagation  constant  when  the 
line  is  lossless,  and 

Xi  =  R'/(kZo)  (5) 

is  the  normalized  resistance/length,  with  =  (L'/C)^/^ 

the  characteristic  impedance  of  the  lossless  line.  The  real 
part  of  y  is  the  attenuation  constant;  the  imaginary  part, 
the  propagation  constant.  Since  the  propagation  constant 
depends  on  R',  finite  resistance  causes  both  attenuation  of 
the  voltage  and  current  waves  as  well  as  phase  shifts. 

The  relationship  between  resistance/length  R'  and  the 
incident  source  power/length  P'=P/L  is  readily  obtained  from 
Eqn  (2)  (Fig.  1) : 
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Figure  2  Model  of  an  Infintessimal 

Section  of  a  Resistive  Line 
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R'  =  Eg/(P'e[iaT) 


(6) 


This  relationship  will  be  used  in  all  subsequent 
calculations . 

In  addition  to  the  obvious  effect  of  reduced  radiation 
efficiency,  finite  resistivity  in  a  long  antenna  changes  the 
radiation  pattern.  First,  if  the  resistivity  is  constant 
along  the  antenna  line,  the  current  decays  from  the  terminal 
towards  the  ends,  effectively  shortening  the  length  of  the 
antenna  and  increasing  the  angular  distribution  of 
radiation.  Secondly,  the  additional  phase  shift  caused  by 
the  resistive  loss  changes  the  phase  relationship  between 
the  radiation  fields  generated  by  the  different  segments  of 
the  antenna,  and  the  radiation  pattern  changes  as  a 
consequence  of  the  changed  interference  pattern. 

Other  antennas  can  be  modeled  as  a  multiple  of  line 
antennas,  so  the  comments  on  resistive  line  antennas  above 
also  apply.  In  particular,  the  sheet  current  in  the  bow  tie 
antenna  studiedd  in  this  work  is  modeled  by  several  line 
loops . 

Boundary  conditions  used  are  zero  current  at  open  ends  of  a 
line,  and  zero  voltage  at  shorted  ends  [2] . 

Finally,  there  is  one  incidental  but  important  difference 
between  a  free-standing  metallic  antenna  and  an  antenna  on  a 
semiconductor  wafer.  Semiconductors  usually  have  very  high 
dielectric  permittivities,  typically  10  times  that  of 
vacuum.  Thus  the  wavelength  in  semiconductors  are  roughly  3 
to  4  times  shorter  than  in  vacuum,  i.e.,  the  physical 
dimensions  of  an  antenna  on  a  semiconductor  are  shorter  than 
those  of  a  free-standing  antenna  of  the  same  electrical 
dimensions  by  a  factor  of  3  or  4.  But  the  radiation  is  into 
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the  vacuum,  so  the  antenna  on  a  semiconductor  has  an 
effectively  reduced  size,  leading  to  lower  efficiency.  By 
dimensional  analysis,  it  can  be  shown  that  the  efficiency  is 
reduced  by  about  the  ratio  of  permittivities,  i.e.,  about  a 
factor  of  10,  which  is  borne  out  by  numerical  calculations. 
The  wave  vector  and  wavelength  in  semiconductor  (vacuum) 
will  be  denoted  as  k  (k^)  and  X  .  The  relationships 
between  these  quantities  are:  k  =  nk^,  X  =  X^/n  where  n  is 
the  refractive  index  of  the  semiconductor.  For  all  the 
calculations  below  for  silicon,  n=3.45  . 

The  general  discussions  above  are  in  terms  of  radiating 
antennas.  By  reciprocity,  the  same  apply  to  receiving 
antennas . 

Formulas  used  in  the  following  calculations  are  collected  in 
Appendix  I.  Calculations  were  carried  out  on  a  personal 
comupter  with  the  software  Mathcad  for  the  linear  antenna, 
and  on  a  Sun  workstation  with  codes  written  in  C  for  the 
longer  procedures  for  the  bow  tie;  C  programs  written  for 
this  contract  are  listed  in  Appendix  II. 

Calculations  are  performed  for  optical  illuminatin  of 
silicon.  Source  power  for  optical  illumination  of 
germanium^  and  electron  excitation  of  silicon  and  germanium, 
can  be  found  by  scaling  using  the  values  in  Table  1. 

(ii)  The  Linear  Dipole 

The  symmetric  linear  dipole  antenna,  illustrated  in  Fig.  3, 
has  been  studied  with  respect  to  its  radiation  pattern, 
resistivity/length,  and  radiation  resistance.  Because  of 
azimuthal  symmetry,  radiation  only  varies  with  0. 
Comparisons  have  been  made  between  a  free-standing  metallic 
antenna  and  one  on  a  semiconductor  wafer. 
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(a)  Lossless  Antennas  For  a  lossless  antenna,  free-standing 
or  on  semiconductor,  the  radiation  pattern,  and  therefore 
directional  gain,  does  not  change  signif icanlty  for  2L  up  to 
one  wavelength  (kL=koL=:7i)  (Fig.  4(a)  and  4(b)).  The 

radiation  resistance,  however,  peaks  at  one  wavelength  (Fig. 
5(a)  and  Fig.  5(b)).  Fig  5  shows  that  the  radiation 
resistance  is  about  10  times  lower  for  the  antenna  on 
silicon  than  free-standing.  For  optimal  performance,  the 
antenna  on  silicon  should  be  made  one  wavelength  long  (1-X), 
which  is  the  simplest  useful  antenna.  Again,  the  wavelength 
is  for  silicon,  therefore  the  physical  length  of  the  1-X 
antenna  is  n  or  3.45  shorter  than  the  wavelength  in  vacuum, 
which  is  still  shorter  than  a  free-standing,  1/2-X  antenna. 
The  following  calculations  will  therefore  be  on  the  l-X 
antenna  on  silicon. 

(b)  Resistive  1-A.  Antenna  on  Silicon,  Uniform  Illumination 

Resistive  loss  is  introduced  in  the  antenna  line  by  Eqn.  (4) 
via  the  dimensionless  parameter  which  is  related  to  the 

resistance/length  and  the  incident  power/length  by  Eqns.  (5) 
and  (6)  .  It  is  assumed  that  the  illumination  on  silicon  is 
uniform,  so  that  the  resistance  is  constant  along  the 
antenna  line.  Fig.  6  shows  the  radiation  resistance  Rj-ad 

the  l->.  antenna  vs  x^.  At  Xi=0.22,  the  has  decreased 
by  a  factor  of  2  from  the  perfectly  conducting  case  (X]^=0), 
i.e.,  the  efficiency  of  the  antenna  has  dropped  by  3  dB.  At 
this  point,  the  total  source  power  needed  to  create  the 
antenna  can  be  found  from  Eqns  (5)  and  (6)  to  be  ~2X2 Watts 
for  optical  illumination  with  X  in  cm.  This  power  level 
agrees  well  with  the  estimate  presented  earlier. 

(c)  Resistive  1-A,  Antenna  on  Silicon,  Non-uniform 
Illumination  The  current  in  the  uniformly  illuminated 
antenna  above  decays  by  1/e  in  a  distance  of  about  l.AX  for 
^1=0. 22,  when  Rrad  drops  by  3  dB,  thus  effectively 
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Figure  5  Radiation  Resistance  vs.  Antenna  Length 
(a)  Free-standing  (b)  Silicon 


shortening  the  antenna  and  reducing  its  efficiency.  The 
obvious  solution  seems  to  be  illuminating  the  end  of  the 
antenna  more  than  near  the  feed  terminal,  i.e.,  higher 
conductivity  as  one  moves  away  from  the  terminal,  more 
uniform  current  distribution  can  be  achieved. 
Alternatively,  for  the  same  R^ad  desired,  less  source  power 

is  needed.  This  is  not  borne  out  by  actual  calculations. 
Non-uniform  illumination  leads  to  a  non-uniform  propagation 
constant,  the  wave  equation  for  which  has  no  known  general 
analytical  solution.  Using  adiabatic  approximation,  valid 
when  the  rate  of  change  of  the  propagation  constant  is  slow, 
calculations  have  been  made  for  (kz) =-0 . 14 (kz) +0 . 44,  whose 

average  kL=7i  is  0.22,  same  as  the  example  above,  and 
therefore  requiring  the  same  total  optical  power.  The 
current  distribution,  radiation  pattern,  and  radiation 
resistance  are  all  almost  indistinguishable  from  the  uniform 
illumination  case  of  X]^=0.22  .  In  retrospect,  the  case  of 

non-uniform  illumination  offers  little  advantage  in 
comparison  to  the  case  of  uniform  illumination  with  the  same 
total  optial  power:  since  the  ends  of  the  antenna  have  to  be 
illuminated  more  to  keep  the  conductance  high,  the  parts 
near  the  terminal  must  be  illuminated  less,  thus  attenuating 
the  current  at  the  very  beginning.  The  only  situation  which 
calls  for  non-uniform  illuminatin  is  when  a  long 
(>wavelength)  antenna  is  desired.  The  source  power  required 
in  such  case  will  be  high. 

(d)  Arrays  of  Linear  Dipoles  Arrays  of  antennas  provide 
higher  gains  than  single  antennas.  The  maximum  directional 
gain  of  a  array  of  N  antennas  is  times  the  maximum 
directional  gain  of  a  single  antenna.  The  maximum 
directional  gain  of  a  single  1-X  antenna  is  only  slightly 
higher  than  unity.  An  array  of  3  provides  increases  the 
directional  gain  by  a  factor  of  9  to  10,  yet  the  optical 
source  power  needed  only  increases  by  a  factor  of  3  to  about 
1  kW,  which  we  consider  practical. 
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(iii)  Bow-Tie  Antenna 


The  bow-tie  antenna,  shown  in  Fig.  7a,  can  be  understood  as 
the  sum  of  many  triangular  wire  antennas  each  with  a 
slightly  different  apex  angle  (Fig.  7b)  .  The  horizontal 
part  of  the  current  contributes  little  to  the  total 
radiation,  thus  each  component  dipole  resonates  at  a 
slightly  different  frequency,  leading  to  the  broad-band 
characteristic  of  the  bow-tie.  In  calculating  the  bow-tie 
antenna  characteristics,  this  model  of  multiple  line  dipoles 
is  used  (Appendix  I).  The  justification  of  this  model  lies 
in  the  good  agreement  between  measurement  [3]  and 
calculation  (Fig.  8).  For  an  antenna  on  semiconductor,  some 
problem  arises  due  to  the  particular  geometry  of  the 
structure.  Uniform  illumination  of  the  antenna  pattern 
leads  to  uniform  conductivity  along  the  current  path. 
However,  the  width  of  the  current  path  increases  along  the 
path,  thus  the  conductance  also  increases,  leading  to  a  non- 
uniform  propagation  constant.  This  problem  exists  even  for 
metallic  antenna:  one  way  to  look  at  it  is  that  as  the  width 
increases,  the  characteristic  impedance  decreases,  again 
leading  to  non-uniform  propagation  constant.  For  this 
study,  we  will  approximate  the  problem  by  using  a  constant 
resistance/length  for  each  of  the  loops  which  consistute  the 
whole  bow-tie.  Justification  lies  on  the  calculated  linear 
dipoles  above  (cf.  Section  (ii-b) ) . 

The  bow  tie-antenna  is  modeled  by  a  series  of  triangular 
wire  antennas  with  the  apex  angle  separated  by  2.5°.  The 
feed  terminal  is  at  the  apex,  and  the  current  is  required  to 
be  zero  in  the  middle  of  the  base  opposite  the  apex  (Fig. 

•  Each  wire  is  assumed  to  have  a  normalized 
resistance/length  x^  (Eqn  5)  .  The  total  radiated  field  is 

calculated  by  adding  the  fields  from  each  wire  vectorially. 
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The  radiation  pattern  and  radiation  resistance  are 
calculated  for  a  metallic  free-standing  bow-tie  to  check 
that  they  agree  with  with  measurements  [3],  which  they  do  if 
a  small  amount  of  loss  is  included,  which  is  not  surprising 
as  resistive  losses  are  present  at  high  frequencies  even  in 
metals.  For  simplicity,  we  will  only  mention,  but  not  show 
any  graph,  that  our  calculations  show  that  the  gain  of  the 
bow-tie  is  comparable  to  that  of  the  half-wave  dipole,  and 
the  radiation  pattern  is  also  similar  except  for  a  slight 
lack  of  azimuthal  symmetry.  The  major  difference  from  the 
half-wave  is  that  as  the  electrical  length  of  the  bow-tie  is 
increased  beyond  a  wavelength,  the  radiation  resistance  will 
not  rise  periodically  to  other  peaks  as  in  the  linear  case. 

The  following  calculations  are  for  a  bow-tie  antenna  with  a 
full  angle  of  60°.  The  radiation  resistance  vs  electrical 
length  kA  is  shown  in  Fig.  9  for  several  values  of 
normalized  resistance/length  X]^.  The  resistance  peaks  near 

the  antenna  length  2A=X  (kA=180°)  ,  as  expected.  Fig.  10 
shows  the  maximum  radiation  resistance  vs  x^.  When  Xi=0.24, 

the  radiation  resistance  has  fallen  by  -3db  from  its 
lossless  case.  To  calculate  the  source  power  needed  to 
induce  this  level  of  resistance,  one  equates  the  total 
resistance  presented  by  the  parallel  wire  antennas  to  the 
current  to  the  total  resistance  presented  to  the  sheet 
current  in  the  bow-tie,  and  the  source  intensity  I  can  be 
shown  to  be  (Fig  7b) : 

I  =  (Eg/e)  (N/2a)ln(A/Zo)/[HT(kA)Zo] 

The  power  is  the  area  of  the  bow-tie  times  I.  For  silicon 
and  kA=7:,  x=0.24,  the  power  needed  to  illuminate  a  circle 
enclosing  the  bow-tie  is  about  1  kW,  again  in  agreement  with 
the  estimate  above. 
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Antenna  Radiation  Resistance  Rr  vs  Length  A 

alpha=30  deq.  n=3.45  for  different  losses  X=R/(KZo) 
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Semi-reconfigurable  Antennas 


Higher  performance  and  lower  source  power  can  be  achieved  by 
trading  off  some  degree  of  reconfigurability  and  stealth. 
Metallic  patterns  can  be  laid  permanently  on  the 
semiconductor,  and  light  or  electrons  will  only  be  used  to 
connect  the  different  parts  of  the  metallic  patterns  to  form 
a  complete  antenna.  For  this  application,  directional 
sources  like  the  laser  and  electron  beams  are  suitable,  and 
the  saving  in  source  power  can  be  considerable,  since  the 
power  required  scales  as  the  square  of  the  distance  between 
2  conducting  patterns  (Eqn.  2) .  However,  the  cross-section 
of  reflecting  an  incident  radar  wave  increases  as  the  area 
of  the  metallic  patterns,  even  if  the  patterns  are  not 
connected. 
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Heating  due  to  Optical  or  Electron  Excitation 


The  heating  due  to  optical  or  electron  excitation  can  be 
estimated  by  considering  the  temperature  rise  after  an 
excitation  pulse.  For  the  applications  considered,  an 
optical  or  electron  pulse  excites  the  semiconductor  to  write 
an  antenna.  Usually  the  pulse  is  much  shorter  than  the 
thermal  time  constants,  therefore  the  excitation  can  be 
regarded  as  an  impulse  whose  energy  deposited  in  the 
semiconductor  causes  a  rise  in  temperature  which 
subsequently  decays  to  the  ambient.  The  maximum  rise  in 
temperature  AT  is  immediately  after  the  exciation  pulse. 
From  the  heat  diffusion  equation  [4], 

AT  =  Q/pC 

where  Q  is  the  energy  density  deposited,  p  the  density  of 
the  semiconductor,  and  C  the  specific  heat.  In  the  worst 
case  of  writing  line  antenna,  assume  a  width  of  5  mm, 
optical  penetration  depth  0.5  pm,  incident  power/length  of 
30  W/cm,  turned  on  for  20  ps,  p=2 . 3  gm/cm^,  C=0.7  J/gm-°C 
(for  silicon),  then  Q=24J/cm^,  and  AT=15°C  which 
subsequently  subsides.  For  electron  excitation,  the  rise  in 
temperature  is  3  times  as  high,  still  insufficient  to  be  of 
concern . 

If  the  excitation  pulse  repetition  rates  increases, 
eventually  heat  cannot  be  conducted  fast  enough.  To 
estimate  the  maximum  rep.  rate,  consider  the  time  heat  takes 
to  diffuse  from  one  face  of  the  silicon  where  it  is  excited 
to  the  opposite  face  at  ambience.  For  a  typical  thickness  s 
of  1/4  mm,  the  diffusion  time  is  s^/D,  where 
D=dif fusivity=0 . 9  cm^/s,  or  0.7  ms,  allowing  a  rep.  rate  of 
1 . 4  kHz . 
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Source  Requirements  and  Available  Sources 


The  only  practical  lasers  are  the  diode  lasers  which  alone 
among  lasers  have  efficiency  of  tens  percent.  Even  they 
have  to  be  cooled  with  thermal  electric  coolers  for  units 
which  put  out  over  100  mW  or  so.  For  pulsed  applications, 
the  energy  needed  is  about  power  x  (silicon  carrier 
lifetime)  or  a  few  mJ.  Commercial,  table-top  diode  laser 
bars  can  deliver  such  energy  levels  at  kHz  repetition  rate. 

For  incoherent  optical  sources,  the  consumer  photographic 
flash  lamp  puts  out  several  joules  in  microseconds.  Even 
accounting  for  losses  in  collimation  and  spectral  mismatch 
with  semiconductors,  the  flash  lamp  delivers  more  than 
enough  energy  for  reconf igurable  antennas  below  1  GHz. 
Repetition  rates  of  photographic  lamp  systems  are  low,  ~  Hz. 
However,  compact  flash  lamp  systems  can  be  constructed  up  to 
kHz  rep  rate. 

The  vacuum  field  emitter  arrays  (FEAs)  can  deliver  kV  and  piA 
per  tip,  with  power  density  in  excess  of  10  kW/cm^  [5]  . 
In  comparison  with  optical  sources,  the  electron  FEA's  are 
more  efficient,  close  to  100%,  and  there  is  little  back- 
scattering  for  the  low  energy  (<  kV)  electrons.  However, 
once  the  electron  enters  the  semiconductor,  only  about  30% 
of  its  energy  is  used  for  generating  charge  carriers.  Thus 
the  overal  efficiency  of  the  two  sources  are  comparable, 
with  the  eleectrons  generating  3  times  more  heat  in  the 
antenna.  Although  reliability  and  costs  of  FEAs  may  be  of 
concern,  the  technology  is  under  rapid  development  and 
driven  by  other  commercial  applications,  which  eventually 
will  make  practical  electron  sources  available  for  this 
application . 
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Thus  the  source  requirements  can  be  met  with  currently 
available  technology,  with  the  cheapest  source,  the 
flashlamp,  most  appealing. 
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Conclusions  and  Recommendations 


Although  detailed  calculations  were  performed  on  only  2 
antenna  structures,  some  general  conclusions  can  be  drawn  as 
more  general  antenna  structures  can  be  viewed  as 
superpositions  of  the  linear  antennas.  First,  because  of 
the  unavoidable  current  decay,  reconf igurable  antennas  of 
dimensions  much  longer  than  a  wavelength  is  not  practical 
unless  for  extremely  high  frequencies.  For  high-gain 
antennas,  arrays  of  smaller  antennas,  each  not  much  larger 
than  a  wavelength,  have  to  be  used  instead.  For  a  single 
antenna  of  a  wavelength  long,  the  lowest  practical  frequency 
is  about  1  GHz,  which  requires  several  watts  of  optical 
power  or  1  kilowatt  of  electron  power.  The  required  power 
scales  as  wavelength  squared. 


Due  to  the  mismatch  of 
semiconductors  and  air,  an 
metallic  or  reconf igurable,  is 
standing  antenna  of  the  same 
actor  of  roughly  the  ratio  of 
for  most  semiconductors. 


dielectric  constants  in 
antenna  on  semiconductor, 
less  efficient  than  a  free- 
electrical  structure,  by  a 
permittivities,  or  about  10 


For  wire-type  antennas,  well-defined  conducting  patterns 
demand  directional  sources  like  lasers  and  electron  beams. 
For  large  aperture  type  antennas,  large  conducting  patterns 
can  use  incoherent  sources  like  the  flash  lamp.  Laser 
systems  of  output  power  over  1  kW  or  so  are  expensive  and 
bulky.  Flash  lamps  of  output  energy  of  several  joules  are 
compact  and  inexpensive,  as  used  in  consumer  photography. 
Electrons  from  vacuum  microtips  are  still  in  the  development 
stage,  and  electrons  do  have  about  only  30%  efficiency  in 
generating  charge  carriers  in  semiconductors. 


For  these  reasons,  it  appears  that  high-performance 
reconf igurable  antennas  are  expensive,  if  indeed  practical 
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at  all,  to  design  and  construct.  The  use  of  reconf igurable 
antennas  should  then  be  confined  to  more  elementary 
structures,  unless  reconfigurability  and  stealth  must  be 
achieved  at  all  cost. 

However,  a  trade-off  can  be  made  between  power  and  the 
degrees  of  reconfigurability  and  stealth  on  the  one  hand, 
and  performance  on  the  other.  If  light  or  electrons  are 
used  only  to  connect  some  permanent,  metallic  patterns  to 
form  antennas,  then  high-performance  can  be  achieved  at 
lower  power  at  a  sacrifice  or  complete  stealth  and 
reconfigurability. 
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Appendix  I  Mathematical  Formulas  and  Equations 

The  following  notations  are  used: 

n  =  index  of  refraction  of  semiconductor 
X,  o=  wavelength  in  vacuum 
JCq  =  wave  vector  in  vacuum 

k  =  nk^  =  wave  vector  in  semiconducto  in  lossless  case 
Zq  =  characteristic  line  impedance  in  lossless  case 

R'  =  resistance/length  of  line 
Xi  =  R'/(kZQ)  =  normalized  resistance 
r\  =  (Pq/Sq)  =  characteristic  impedance  of  vacuum 
=  376  Q 

The  following  formulas  and  equations  can  be  found  in  [2], 
adapted  for  this  particular  study.  The  propagation  constant 
Y  is  defined  in  Eqn  (4) : 

y  =  ( jcoC'R'  -  co2l'C'  )  1/2 

=  jk(l  +  x^)  i/2exp  (-0 . 5jtan~ix2)  (Al) 

Now  consider  a  current  I(z')  flowing  in  a  straight  line 
along  z'  which  lies  on  the  x-z  plane  with  an  angle  b  from 
the  z-axis  (Fig.  A-1) .  The  radiation  fields  can  be  derived 
from  the  vector  quantity  N[2]: 

N  =  [z  cos(b)  +  X  sin(b)]  J  I (z ’ ) exp ( j k^z ' cos (p) )  dz ' 

{A2) 

where 

cos  (p)  =  z'  r 
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For  linear  dipoles  of  length  L  on  the  z-axis,  b-0,  p-6, 

boundary  conditions  I(z'=+L)=0.  Then  for  input  current  of 
amplitude  Iq,  the  current  distribution  is: 

I^[e-YZ'  -  eYz'-2yL]  z'>0 

I(z')  = 

I^feYZ'  _  e_Yz'-2YL]  z'<0 

With  the  normalized  integration  variable  s: 
s  =  kz' 

and  normalized  propagation  constant  g(s): 
g(s)  =  yz ' 

N  can  be  written  as 

kL 

N  =  (zl^/k)  [  jo[e-Y2*  _  eyz’-2YL]ejscos(e)/n  ds  + 

0 

j_kL[e"l'2'  _  eYz'-2YL]e:)s'=°s(®)/^  ds  ] 


=  (zlo/k)  F(kL,n,0)  (A6) 

which  defines  the  expression  F.  The  total  power  radiated 
into  all  solid  angles  is 

W  =  ri  Io2/ (2  (47in)  2)  {4^  dO  sin20  [F  (kL,  n,  0]  2 

=  Io2/(2Rrad) 

which  defines  the  radiative  resistance  Rrad- 
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The  bow-tie  antenna  is  modeled  by  a  number  of  loop  antennas 
with  b's  separated  by  2.5°(Fig  A2) .  The  boundary  conditions 
are  that  the  currents  vanish  in  the  middle  of  the  horizontal 
(parallel  to  x-axis)  segments.  For  each  straight  segment  of 
each  loop,  the  angle  b  changes  from  b  to  90°  (horizontal 
part)  to  -b.  The  function  F  is  calculated  by  summing  all 
the  F's  for  each  loop.  The  integral  in  Eq  (A7)  now  depends 
on  (j)  as  well. 
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Appendix  II  Program  Listing  of  Bow-Tie  Calculations 

/*  bow_tie  antenna  plate  simulation  */ 

/*  radiation  resistance  parameters  :  x,  y,  alfa  */ 

#include  <stdio.h> 
iinclude  <math.h> 

#define  Nz  25 
#define  N  200 
#define  n  3.45 
#define  pi  3.141592654 
idefine  dALFA  2.5 
#define  ALFA  30. 

/*  -o  c05.out  */ 
main ( ) 

{ 

FILE  *fopen( ) , *pout; 

double  theta,dtheta,Kr,Rr, fi,dfi; 

double  Nx( ) ,Nxi ( ) ,Nzz ( ) ,Nzi ( ) , N3r ( ) , N3i ( ) ; 

double  suml, sum2 , si , s2 ; 

double  alfa [37] , y [37 ] , Ntheta [ 37 ] , Nf i [37 ] , Nthetai [ 37 ] ,Nfii [37] ; 

double  cl, c2 , c3 , c4 ; 

double  A, A1 , wavelength, xl ; 

double  Am,Rm,M; 

int  j 1 , j2 , j , q, NO , i ; 

M=  ALFA/dALFA+  0.5; 

M=f loor (M) ; 

Am=0 . ;  Rm=0 . ; 

wavelength=0 . 5 ; / *f =600MHz*/ 
for  { i=l ; i< (M+1 ) ; i++ ) 

{ 

alfa[i] =dALFA*i*pi/360; 

) 

pout  =  f open  (  '' nx3  0R3  0 OAl .  dat " ,  " w" )  ; 
xl=0 .24; 

A=15.  ; 

for  (q=l;q<=l;  q++) 

{ 

Al=A*wavelength/360 ; 
for  { i=l ; i< (M+1 ) ; i++ ) 

{ 

y [i] =A1/ (wave length*cos (alfa [i] ) /4) ; 

} 

N0=N/2; 
dtheta=pi/N0 ; 
dfi=2*pi/N; 
suml=0 . 0 ; theta=0 . 0 ; 
for  ( jl=0; jl<=(N0/2) ; jl++) 

{ 

if((jl%2)  ==0)sl=2.0; 
else 
s  1  =  4  .  ; 

if(jl==0  II  jl==N0/2)sl=l. ; 
sum2=0 . 0 ; f i=0 . 0 ; 
for  (j2=0; j2<=(N/4) ; j2++) 

{ 

if((j2%2)  ==0)s2=2.0; 
else 
s2=4. ; 

if(j2==0  II  j2==N/4)s2=l. ; 
for  ( i=l ; i< (M+1 ) ; i++) 

{ 

Ntheta[i]  =  (Nx(xl,y [i] , alfa [ i] , theta, f i ) +N3r (xl , y [ i ] ,alfa[i]  ,  theta,  fi) ) 
Ntheta [ i ] =Ntheta [ i ] *cos (f i) *cos (theta) ; 

Ntheta [i]=Ntheta[i]-Nzz(xl,y[i] ,alfa[i] , theta, fi) *sin (theta) ; 

Nfi [i] =(Nx(xl,y [i] , alfa [ i] , theta, f i) +N3r (xl,y [ i] , alfa [ i] , theta, fi) ) *sin(fi) 
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Nthetai [i]=(Nxi(xl,y[i] ,alfa[i] , theta, f i ) +N3 i(xl,y[i] ,alfa[i] , theta, f i ) ) ; 
Nthetai [ i ] =Nthetai [ i ] *cos ( f i ) *cos (theta) ; 

Nthetai [ i ] =Nthetai [ i ] -Nzi (xl ,y [ i ] , alfa [ i ] , theta, f i ) *sin ( theta) ; 

Nf ii [ i ] = ( -Nxi (xl , y [ i ] , alfa [ i ] , theta, f i ) +N3i (xl,y[i] ,alfa[i] , theta, f i ) ) *sin ( f i ) ; 

} 

cl=0 . ;  c2=0.;  c3=0.;  c4=0.; 
for  ( i=l ; i< (M+1 ) ; i++ ) 

{ 

cl=cl+Ntheta [ i ] ; 
c2=c2+Nthetai  [ i ]  ,- 
c3  =c3+Nf i [ i ] ; 
c4=c4+Nfii [i] ; 

} 

Kr= { cl*cl+c2 *c2+c3 *c3+c4*c4 ) ; 
sum2=sum2+s2*Kr*sin (theta) ; 
f i=f i+df i ; 

} 

sum2=sum2*df i/3 ; 

suinl  =  suml  +  sl*sum2  ; 
theta=theta+dtheta ; 

) 

Rr=30*pi*suml*dtheta/3 ; 

Rr=Rr/ ( (M) * (M) ) ; 

Rr=Rr/2; 

Rr=Rr*8; 

fprintf  (pout,  "%e  %e\n'',A,  Rr)  ; 
printf("%e  %e\n",A,  Rr) ; 

A=A+5 . ; 

} 


double  Nx (xl , y, alfa, theta, fi) 
double  xl , y, alfa, theta, fi ; 

{ 

double  X, z , dz , V, u , sum, s,J,Jl,pl,nl, f; 

int  i  ; 

z=0.1; 

sum=0 ; 

dz=l . /Nz ; 

for  ( i=0 ; i<=Nz ; i++ ) 

C 

x=xl; 

f =atan (x) /2 ; 

v=n*0 . 5 *pi *pow ( ( l+x*x) , 0 . 25 ) *sin ( f ) ; / *  v=0.5*pi*a  */ 

u=n*0 . 5 *pi *pow ( { l+x*x) , 0 . 25 ) *cos ( f ) ; / *  u=0.5*pi*b  */ 

nl  =  0 . 5*pi*y*sin (alfa) *sin (theta) *cos (f i)  ; 

pl=0 . 5*pi*y*cos (alfa) *cos (theta) ; 

if((i%2)  ==0)s=2.0; 

else 

s  =  4  .  ; 

if(i  =  =  0  I  I  i  =  =Nz)s  =  l.; 

J=cos (u*y* (1+sin (alfa) ) ) *cos (u*y*(l+sin(alfa)-(z)))*sinh(v*y*(l+sin(alfa)-(z))); 
J=J+sin(u*y*(l+sin(alfa) ) )*sin(u*y*(l+sin{alfa)-(z) ) ) *cosh ( v*y* (l+sin(alfa)-(z))); 
sum=sum+s*y*exp (-v*y* {l+sin{alfa) ) ) *J*sin(nl* (z) ) *sin(pl* (z) ) ; 
z=z+dz ; 

} 

Jl  =  -sum*dz/ (3 . )  ; 
return (J1 ) ; 

} 

double  Nxi (xl , y, alfa , theta, fi ) 
double  xl , y, alfa , theta, fi ; 

{ 

double  X,  z ,  dz ,  V,  u,  sum,  s ,  L,  LI,  pi,  nl ,  f  ,- 
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int  i  ; 
z=0 . 1 ; 
sum=0 ; 
dz=l . /Nz ; 

for  ( i=0 ; i<=Nz ; i++) 

{ 

x=xl  ; 

f =atan (x) /2; 

v=n*0 . 5*pi*pow( (l+x*x) , 0 .25) *sin (f ) ; /*  v=0.5*pi*a  */ 

u=n*0 . 5*pi*pow( (l+x*x) , 0 .25) *cos (f ) ;/*  u=0.5*pi*b  */ 

nl=0.5*pi*y*sin(alfa) *s in (theta) *cos (f i) ; 

pl=0 . 5*pi*Y*cos (alfa) *cos (theta) ; 

if((i%2)  ==0)s=2.0; 

else 

s=4 .  ; 

if(i==0  i|  i==Nz)s=l.; 

L=cos (u*y* (l+sin(alfa) ) ) *sin (u*y* ( 1+sin (alfa) - (z) ) ) *cosh(v*y* (1+sin (alfa) - (z) ) ) 
L=L-sin(u*y* (l+sin(alfa) ) ) *cos (u*y* ( 1+sin (alfa) -z) ) *sinh (v*y* ( 1+sin (alfa) -z)); 
suin=sum+s*y*exp(-v*y* (1+sin (alfa) ) ) *L*sin (nl* ( z) ) *sin(pl* (z) ) ; 
z=z+dz ; 

} 

Ll  =  -sum*dz/ (3 . )  ; 
return (LI ) ; 

) 


double  Nzz (xl,y, alfa, theta,  fi) 
double  xl , y, alfa, theta, fi ; 

( 

double  X, z, dz, V, u, sum, s, J, Jl,pl, nl, f ; 

int  i  ; 

z=0 . 1 ; 

sum=0 ; 

dz=l . /Nz ; 

for  ( i=0 ; i<=Nz ; i++) 

{ 

x=xl  ; 

f=atan (x) /2 ; 

v=n*0 . 5*pi*pow( ( l+x*x) , 0 .25) *sin (f ) ; /*  v=0.5*pi*a  */ 
u=n*0 . 5*pi*pow( (l+x*x) , 0 .25) *cos (f ) ; /*  u=0.5*pi*b  */ 

nl=0 . 5*pi*y*sin (alfa) *s in (theta) *cos (f i) ; 
pl=0 . 5*pi*y*cos (alfa) *cos (theta) ; 
lf((i%2)  ==0)s=2.0; 
else 
s  =  4  .  ; 

if(i==0  II  i==Nz)s=l.; 

J=cos (u*y* (1+sin (alfa) ) ) *cos (u*y* ( 1+sin (alfa) -( z ) ) ) *sinh (v*y* ( 1+sin (alfa ) - (z) ) ) 
J=J+sin(u*y* (1+sin (alfa) ) ) *sin (u*y* ( 1+sin (alfa) - (z) ) ) *cosh (v*y* ( 1+sin (alfa) - (z) ) ) 
sum=sum+s*y*exp(-v*y* (1+sin (alfa) ) ) *J*cos (nl* (z) ) *cos (pi* (z) ) ; 
z=z+dz ; 

) 

Jl=sum*dz/ (3 . ) ; 
return ( J1 )  ; 

} 


double  Nzi (xl,y, alfa, theta, fi) 
double  xl,y, alfa, theta, fi; 

{ 

double  X, z, dz, v,u, sum, s , L, LI , pi , nl , f ; 
int  i  ; 

Z=0.1; 
sum=0 ; 
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dz=l . /Nz ; 

for  ( i=0 ; i<=Nz ; i++) 


{ 

x=xl  ; 

f=atan (x) /2 ;  .  «  c*  ; 

v=n*0 . 5*pi*pow( (l+x*x) , 0 .25) *sin(f ) ; /  v-0.5  Pi  ^ 

u=n*0 . 5*pi*pow( (l+x*x) ,0.25) *cos (f ) ; /*  u=0 . 5*pi  b 

nl=0.5*pi*y*sin(alfa) * sin (theta) *cos ( f i) ; 
pl  =  0 . 5*pi*y*cos (alfa) *cos (theta) ; 
if((i%2)  ==0)s=2.0; 

else 


*/ 

*/ 


sum=sum+s*y*exp(-v*yMl+sin(alfa)))*L*cos(nlMz))*cos(pl  (z)), 

z=z+dz ; 


) 

Ll=sum*dz/ (3 . ) ; 
return (LI ) ; 

] 


double  N3r(xl,y, alfa, theta, fi) 
double  xl , y , alfa , theta , fi ; 

double  X, z,dz, v,u, sum, s , L, LI , pi , nl , f , zl ; 

int  i  ; 

z=0 . 1 ; 

sum=0 ; 

dz=l . /Nz; 

for  (i=0;i<=Nz;i++) 


( 

x=xl  ; 

f=atan (x) /2 ;  ^  c*  ■* 

v=n*0.5*pi*pow((l+x*x),0.25)*sin(f);/*  v=0.5*pi  a 

u=n*0 . 5*pi*pow( (l+x*x) , 0 .25) *cos (f ) ; /*  u-0 . 5  pi  b 

nl=0.5*pi*y*sin{alfa)*sin (theta) *cos ( f i ) ; 
pl=0 . 5*pi*y*cos (alfa) *cos (theta) ; 
if((i%2)  ==0)s=2.0; 

else 


*/ 

*/ 


s  =  4  .  ; 

if(i==0  II  i==Nz)s=l.; 

zl=  (1+sin  (alfa)  *  (1-z)  )  ;  •  , 

L=cos(u*y* (l+sin(alfa) ) ) *sin (u*y * ( 1+sin (alfa) - ( zl ) ) ) , 

T.-L*cosh  ( v^v*  ( 1-*- s  in  ( a.1  f  a )  -  ( zl )  )  )  ;  . 

L=L-.ln(u:y-(l*.in(L£.||)-co.Wy(l«in(.lfa>-»l))*s;nMVy(l*  =  in(alEa)- 

Ll=s*v*exp ( -v*y* ( 1+sin (alfa) ))*( 0 . 5  sin{alfa))  L, 

Ll=Ll*sin ( 0 . 5*pi*y*sin (theta) *sin(fi)*sin(alfa)*zl) ; 


sum=sum+Ll ; 

) 

Ll  =  -sum*dz/ ( 3  .  )  ; 
return (LI) ; 

] 

double  N3i(xl,y, alfa, theta, fi) 
double  xl,y, alfa, theta, fi; 

double  X, z, dz, V, u, sum, s , J, J1 , pi , nl , f , zl ; 

int  i  ; 

z=0 .1; 

sum=0 ; 

dz=l./Nz; 

for  ( i  =  0 ; i<=Nz ; i  +  + ) 

{ 

x=xl  ; 


-(z)  )  ) 
- (z)  )  ) 


zl)  )  ; 
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f =atan (x)/2; 

v=n*0.5*pi*pow( (l+x*x) , 0.25) *sin(f ) ;/*  v=0.5*pi*a  */ 

u=n*0 . 5*pi*pow ( ( l+x*x) , 0 . 25) *cos (f ) ; /*  u=0.5*pi*b  */ 

nl=0 . 5*pi*Y*sin (alfa) *sin (theta) *cos ( f i) ; 

pl=0 . 5*pi*Y*cos (alfa) *cos (theta) ; 

if((i%2)  ==0)s=2.0; 

else 

s  =  4  .  ; 

if(i==0  II  i==Nz)s=l.; 
zl=(l+sin(alfa)*(l-z) ) ; 

J=cos (u*y* {l+sin(alfa) ) ) *cos (u*y* (l+sin(alfa) - (zl) ) ) *sinh(v*Y* ( 1+sin (alfa) - ( zl ) 

)  )  '• 

J=J+sin (u*y*(l+sin(alfa) ) ) *sin(u*y*(l+sin(alfa)-(zl) ) ) *cosh (v*y* (l+sin(alfa)-(zl) 

)  )  ; 

suni=sum+s*y*exp(-v*Y* (1+sin (alfa) ) ) * ( -0 . 5*sin (alfa) ) *J*sin ( 0 . 5*pi*y*sin ( theta ) *sin( 
f i ) *sin (alfa) *zl ) ; 

z=z+dz ; 

} 

Jl  =  -sum*dz/ (3 . )  ; 
return ( Jl) ; 
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